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Abstract 

Tree amplitudes of the production of two kinds of scalar particles at threshold from 
one virtual particle are calculated in a model of two scalar fields with 0(2) symmetric 
quartic interaction and unequal masses. These amplitudes exhibit interesting factorial 
and exponential behaviour at large multiplicities. As a by-product we observe that 
the kinematically allowed decay of one real particle into n real particles of another 
kind, all at rest, has zero tree amplitude in this model for n > 2. 







1 Introduction 



Recently, considerable attention have been paid to the study of multiparticle pro- 
duction processes both in the instanton sector [|l], ^, ^ and in the conventional 
perturbation theory p, ^ ^. In the latter case, most of the results has been ob- 
tained in simple theories of one scalar field (or p scalar fields with 0{p) symmetry), 
where a number of exact expressions have been derived, e.g., for the amplitudes of the 
production of n particles at the multiparticle threshold by one virtual particle, at the 
tree level |^, ^ and in the one loop approximation |TD|, |TI], |T^. These amplitudes 
grow like n\ at large n, and this non-trivial property is expected to persist to all orders 
of the perturbation theory |T^, |T^, p!2| . 

For calculating the tree amplitudes at the threshold, two methods have been em- 
ployed. One of them is based on recursion relations between diagrams with different 
numbers of final particles 0], and the other makes use of classical field equations 
with special boundary conditions @] . Both of these methods have been generalized to 
incorporate loop corrections [Tl|, |l^ and both work nicely in the simplest theories. 

In this paper we study the tree amplitudes at the multiparticle threshold in the 
theory of two real scalar fields with the lagrangian 

L = lid.^r + lid.x? - ^v^^ - - A(^^ + xr (1) 

The interaction term is 0(2)-invariant, while the mass terms explicitly break 0(2) 
symmetry at mi ^ 1112. This model may be viewed as a testing ground for different 
calculational techniques, and also has some features absent in the simplest theory 
of one scalar field. For example, the model contains an additional dimensionless 
parameter, mi/m2, and one can study the dependence of the amplitudes on this 
parameter. Also, the class of the tree amplitudes is wider than in the theory of one 
field because two kinds of particles can be emitted. Finally, in the case of broken 
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reflection symmetry, ip —(f, there emerges a kinematically allowed possibility (at 
> 1712 > "^i) of ^ decay of an on-s/ie// <y9-particle at rest into x-particles, all of which 
are also at rest (of course, an interesting quantity is the amplitude of this decay and 
not the decay width which is zero because of empty phase space). 

We find that the classical solution method is convenient for the evaluation of 
the tree amplitudes at the threshold in the model (|T|). Due to a special symmetry 



of the corresponding classical equations |]T5[, the explicit solution of the relevant 
boundary value problem can be found, and the explicit form of the tree amplitudes 
at particle threshold is obtained in this paper for the processes riiip + and 
X — ^ Uiif + n2X) where the initial particle is off shell, for both unbroken (< ip >= 0, 

< X >= 0) and broken (< ip >^ 0, < x >= 0) cases. Not surprisingly, these 
amplitudes grow factorially at large n\ and n2, but besides the expected factorials, 
the amplitudes contain exponential terms depending non-trivially on m^jm^ and 
ni/n2- Less obvious is the fact that the tree amplitudes of the decay of an on-shell 
(/9-particle at rest into n2 x-particles with zero momenta (kinematically allowed in the 
spontaneously broken case) is zero for n2> 2. 

In Sect. 2 we obtain the general classical solution to the classical field equations for 
space-independent fields. In Sect. 3 we evaluate the tree amplitudes at the threshold 
in the unbroken symmetry case < (p >=< x >= 0- The broken symmetry case 

< ip >7^ 0, < X >= is considered in Sect. 4. Sect. 5 contains concluding remarks. 



2 Solution to classical equations 

We consider the tree amplitudes of, say, the processes (p — >• ni(p+n2X, where the initial 
particle is off-shell and the final particles are on-shell and have zero spatial momenta. 
According to ref.[|l, the generating function for these amplitudes, ipc{zi,Z2), is the 
solution to the classical field equations for space-independent fields, obeying the con- 
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dition that 

if-^ zi + ifo, X-" Z2 (2) 

as A — > 0, where 

= aie^'"^*, Z2 = aae*"^^* (3) 

and are the particle masses, ipo =< ip > (hereafter we assume < x >= 0), 
and ai and 02 are arbitrary constants. The amphtudes are then 

K,,n, = ^ ^ Vc{z,, ^2) ^^^^^^^ (4) 

Thus, we have to consider the hamiltonian for space-independent fields (we set 
the spatial volume equal to 1; the dependence on the volume can be easily restored), 

H = li^? + lix? + + ^X' + A(v.^ + xr 



It has been found in ref. [0 that this hamiltonian possesses a non-trivial symmetry 
and thus is integrable. It has been also pointed out ||16[ that the system is separable 
in elliptic coordinates. So, we introduce new variables ^{t) and i]{t), 



X = (^^Vj ^ = const. 

By choosing 



a 



2 mi — 1712 



2A 

we find that the variables indeed separate in the Hamilton- Jacobi equation and we 
find for the truncated action S{C,,ri), 



ie - i)(f )' + 2a^{xar +Cy- 2Xaye - Cy^' - + E)e} 

iv' - l){^y + 2a^{AaV +Cy- 2Aa2)aV - Cy^' - Aa^ + E^} 



2(tV (5) 
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where E is the classical energy and J is another integral of motion. Therefore, the 
truncated action is 

S = F{0 + F{r]), 

where 



2a^J - 2a2^2[^^4^4 + (^2/2 _ 2X(x^)a^^^ - (m?aV2 - Xa^ + E)] 



The action function is then 

Stot = -Et + F{0 + F{7]). 

So, the general solution is determined by 

dStot _ , , dFjO , dFjr^) _ ^ 
dE - ^ dE ^ dE ~ °' 

dStot _ n 

where P and to are arbitrary constants. Thus, the general solution is finally obtained 
in the implicit form 

C2 2 

+ [dv^ = t-to, (6) 

■J giO m) 
1 o /■. 1 



where 



/^^^ + ^' /rf7^-- = /3, (7) 



1 1 
mi „ , „\ „ /mi 



Since we are interested in complex solutions, the arbitrary constants E, J, (3 and to 
are also, in general, complex. We will see, however, that in both broken and unbro- 
ken symmetry cases, the relevant solutions correspond to J = 0, E = H{< ip >) = 
(vacuum energy), while to and P are determined by the conditions (13). 



3 Unbroken reflection symmetry 

Let us first consider the case mf > 0, > when < ip >= ipQ = 0. To ensure the 
conditions (Q) we note that in the hmit A — > 0, the fields ip and x do not interact with 
each other, and their energies, and E^, are the two integrals of motion. So, the 
integrals of motion of the non-linear problem, E and J, may be expressed through 
E^p and E^ in this limit. Obviously, E = E^ + E^, and it is straightforward to see 
from eq.(^) that J = —E^ at A ^ 0. 

For positive frequency solutions, eqs. (0) and (|^), one has E^ = E^ = 0, so the 
relevant solution to the non-linear equations is determined by eqs. (|^), (|^) and 
with 

E = J = 

Then the system (|), (|^) can be solved explicitly, and the solution obeying eqs. @ 
and is (see Appendix A for details) 

where 

mi — 7712 

K = 

mi -|- m2 

Notice that at ^2 = 0, the function (P) coincides with the known solution for one 
scalar field P], as it should. 

The amplitudes ip riiip + n2X are given by eq. (^. They are non-zero at 

ni = 2p+ 1, 

n2 = 2q, 
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where p and q are integer. We find 



mTr.^^ ' k) ^''^'^ V.ip-l)liq-l)l ^'^'-^'-^''^^'^ 

(11) 

The finite sum in eq.(|l^) can be transformed into the Jacobi polynomials of the 
argument (1 — 2k^) (see Appendix B), so that the expression for the amplitude takes 
the form 



_ (2p + l)K2g)!.Axp+, ,,, P + g-l ^^ 
~ m'^ml' ^2) ^ ^^^^^ 

(p + g _ P(V-P'0)(1 - 2k') - (g - 1)(1 - «:)«:2PS!/-^'^)(1 - 2k')\ (12) 



An equivalent form of the amplitude can be obtained by transforming the sum in 
eq.(|Tl|) into the hypergeometric function F (see Appendix B). At g > p we have 

_ (2p+2)K2g)!/V+, 

^P>9 2p 2q [r, ) ^ 

17111712 

{/.^^(-ir j^~_^^'^J pF(-p,g;g-p (13) 

while the expression at g < p is obtained by interchanging p <-> g in the curly brackets 
in eq. ([T3|) . 

Eq.(|Tl|) (or, equivalently, eqs. (p!2|) or (0)) is the exact formula for the tree ampli- 
tudes. Its asymptotics at large p and g and p/q = fixed has the form 

(2p+l)!(2g)! ,^,rni p.^ 

Ap,, '-J^^^l^ exp[(p + g)G(— ; q] (14) 

nil 1712 "^2 Q 

The explicit expression, given in Appendix C, is not very illuminating. Notice that 
Ap q grows factorially as could have been expected. 

A particularly simple asymptotics emerges at p ^ cxd, g —>■ oo, p/q 1. In that 
case we obtain (see Appendix C) 



nil m2 



where C{p) depends also on ini/ni2 and grows like some power of p. 



4 Broken reflection symmetry 

At TTif < 0, 777,2 ^ symmetry (p — > —if is spontaneously broken, and 

Imil 

The masses of excitations around this vacuum are 



rriu 



\/2\mi\, = J \mi\'^ + 7712 (16) 



The conditions (0) and (|) are satisfied when the integrals of motion are equal to 

J = 0, E = ——rr = Eyac 

loA 

The explicit solution to the equations of motion, that obeys the conditions (|]) and 
(^, is obtained in the same way as in Sect. 3. We find 



/ Zi 2A 2 2777^ — 2\ 



/i ^1 2A ^2 I ^ 2m^-m^ ^ ^2\-i ..7^ 



/ 2777^ - Zi X 

X = Z2[l - ( ^ j^)^ 



^ 2(y9o 4777^ -m'^ ^ ipo (2777^ + 777,^)3 

The field (|l^ reduces to the known solution for spontaneously broken theory 
of one scalar field [§] at Z2 = 0. On the other hand, eq.(|l8D does not coincide at 
Zi = with the solution for the theory of one field with unbroken symmetry [^]; in 
the diagrammatic language this corresponds to the existence of the diagrams of fig. lb 
absent in the theory of one field. 

In the case of broken symmetry we should distinguish between the processes (f — > 
ni(f + 772X and x + f^2X- The amplitudes of the former, which we denote by 



^ni,2q (the number of x-particles in the final state is even), are generated by (p{zi,Z2). 
The amphtudes of the latter, A^^^2q+i: generated by x{zi,Z2)- We find for the 
process ip riiip + 2qx, 

where 

^ _ 2m^ - 
2m^ + 

In analogy to Sect. 3, the amplitudes in eq. (p!9[) can be expressed through the 
Jacobi polynomials or hypergeometric function. We present here the latter form only, 
which reads 

nr\{2q)\ ^ ^^^^^^^ 

{(-irC^"^ _ f^),^^, g + + l)-^F(-ni + l, g; g-^ + l; 1)]} 

(20) 

Here it is assumed that g > ni; at g < rii the expression for the amplitude is obtained 
by substituting g ^ ni in curly brackets. 

From the latter expression and eg. (PPP , the asymptotics of the amplitude at 
rii ^ oo, g ^ oo, ni/g ^ 1 can be obtained, 

« ...-^f'"" ,VTx)-*"-^[ ^^^\^^Y Cin.,,) (21) 

where C{ni, g) has power-like behaviour in this limit. The asymptotics in the general 
regime, ni —>■ oo, g oo, rii/q =fixed, can also be found and has the general form of 
eq.dH). 

From eq.(0) we find the amplitudes of the processes x ~^ ^if + (2g + 1)X) 
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This expression can again be rewritten in terms of the Jacobi polynomials or hyperge- 
ometric function. At rii oo, q oo, rii/q ^ 1 we have the following asymptotics, 
ni!(2g + l)! /— r(v/2^ + yrrv)^^ 2ni 



We observe that both A^_^ 2q+i ^ni,2q increase factorially at large number of final 
particles. 

To conclude this section, we point out that at certain values of |mi/m2| such that 
f^tfi = n2m^ = 2qm^, there exists a kinematically allowed possibility of the decay of 
a (/9-particle at rest into n2 = 2q ^-particles, all of which are also at rest (according 
to eq.(0), this requires > ml > mf). Eq.([T9|) tells us that this amplitude is in fact 
zero at the tree level, unless n2 = 2 (g = 1). Indeed, the amplitude ^0,2? corresponds 
to the diagrams like those shown in fig.2, where the propagator of the initial particle 
is included. So, the decay amplitude A'^f^2qx obtained from ^0,2^ by truncating the 
leg corresponding to the initial particle, 

K^%x = (VmJ - ml)Al2q (23) 

At g = 1, the amplitude of eq.([T9|) has a pole at = 4g^m^ = 4m^, and one recovers 
the obvious result. 



On the other hand, at g > 1 (^2 > 2), the amplitude of eq.(^) does not have a pole 
at = 4g^m^, and the decay amplitude ( p3D vanishes. This means, for example, 
that the diagrams of fig.2 with truncated initial legs, cancel each other for on-shell 



initial and final particles, all at rest (for the decay (f 4x this fact can of course be 
checked by the direct calculation of the diagrams). 



5 Conclusion 

The model of two scalar fields with 0(2) symmetric interaction and unequal masses 
has the property of integrability and separability of its classical equations for spatially 
homogeneous fields. This property enabled us to obtain the explicit formulas for 
the tree amplitudes at particle threshold for both unbroken and broken reflection 
symmetry (f —ip. These amplitudes are expressed through the Jacobi polynomials 
or, equivalently, finite hypergeometric series. 

The behaviour of the tree amplitudes of the creation of rii t^^-particles and n2 x- 
particles at threshold at large rii, n2 and fixed n\jn2 is suggestive. It has the factorial 
growth and is exponential in the free parameters m^p|m^ and nxjn^^ 

An„n, oc n,\ A(-i+"^)/2 exp[(ni + n2)G(^, ^)] (24) 

up to well defined powers of mi and m2 and pre-exponential factors growing power- 
like at large ni, n2- The function G{^, ^) is calculable at the tree level; its explicit 
form at arbitrary ni/n2 in the unbroken case can be read off from the eq. (^), while 
at ni/n2 — ^ 1 the corresponding expressions are particularly simple, see eqs. (pISD, 
(^ and (^). 



Eq. (^) indicates that there may exist a semiclassical-type procedure for cal- 
culating the leading exponential behaviour of the amplitudes at large rii, n2. Such 
a procedure is explicit in the recent study [1^ of the amplitudes in the theory of 



one scalar field with broken symmetry. So, one may anticipate that this feature is 
characteristic to a wide class of bosonic theories. 

A peculiar property of the model with < ip >^ is that the decay of a (/^-particle 
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into 2q ^-particles at rest, although kinematically allowed at appropriate values of 
particle masses, has zero tree amplitude for 2q > 2. This fact may or may not be 
related to the integrability of the classical theory of spatially homogeneous fields. We 
hope to discuss the cancellation of the corresponding diagrams (that occurs also for 
some other processes) in future. 

The authors are indebted to D.T. Son, E.E. Tareyeva and P.G. Tinyakov for 
helpful discussions. One of us (V.R.) thanks M.B. Voloshin and L. McLerran for 
stimulating conversations at the initial stage of this work. The work of M.L. and S.T. 
is supported in part by the Weingart Foundation through a cooperative agreement 
with the Department of Physics at UCLA. 



Appendix A. Classical solution 

The integrals in eqs.(P), (|^) are straightforward to evaluate at E = J = 0. The 

system (D), (0) then takes the form 

(■0 + mi)(u + mi) r / XT 

)^ '-^ ^ = exp[T2zmi{t - to)] 

[ip - mi)(/i - mi) 

exp[^2im2(t - to - 2/3(t^)] (25) 



(■?/' + m2)(/i + m2) . , 2 



{if) - m2)(/i - m2) 



where if) = \J a^'^ + rn^^ /i = \J aff + m2, a = ml — ml- We note that the original 
fields ip and x can be expressed in terms of /i and ip, 

= ^(^^ - rnl){f/ - ml) 

Therefore, to find the fields (p and x one does not need to evaluate ^(t) and rjit), but 
one just has to solve the algebraic equations (^3]). This is easily done in terms of 
variables (V^/i) and (■?/' + /i). We find that the correct behaviour (||),(0) corresponds 
to the upper sign in eq. (|25|), and that the solution is given by eqs. (^, ([T0|). 
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Appendix B. Reduction to Jacobi polynomials and 



hypergeometric functions 

The sum in eq.(|lT]) can be rewritten in the following form, 



E = E J,, [(P + ? - - (1 - «^)^] TT X 
I y- y- '■• 

[(-p+l)---(-p + 0][(-g + l)---(-g + /)] 
[-{p + q-l) + l]---[-{p + q-l) + l] 

This should be compared with the hypergeometric series |T^, |18[ 

^ ' ' ' ^ " r(a)r(6) ^ r(c + o i\ 

The direct comparison gives 

E = ^^^f 7^^' [(p + g - Kq)F{-p + 1, -g + 1; -p - g + 2; 
q\ 

(1 - «:)«:^ t^i^)izi±l)F(_p + 2, -g + 2; -p - g + 3; 
{-p-q + 2) 



Eq. (0) is then obtained by making use of the relations [T8| between the Jacobi 
polynomials and hypergeometric functions. 

The representation (|l^) is obtained by changing the variable in the sum, / p — l, 
and then by using essentially the same trick. 



Appendix C. Asymptotics of the amplitudes 

Let us first evaluate the asymptotics of the Jacobi poljTiomials P^°''^\x) at r ^ cxd, 
a — oo, a/r = —6 = const. We make use of the generating function for the Jacobi 
polynomials, 

jr(^) = 2°+^i?-^(l -y + + y + R)-^ 
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where R = \/\ — Ixy + y^, \y\ < 1. The Jacobi polynomials are 



XI 



= -Jfe-yyHy)dydy 



(26) 



The integral in this equation is evaluated in the saddle point approximation. The 
saddle point is 



r 

y 



In this way we obtain the asymptotics 



(x + l)(x + 1 + 



(0 - 2y 



p(»'f^\x) = 2V(1 -y + RT Cir, a) 



(27) 



(2^ 



where y and R{y) are taken at the saddle point (|27|) , and C{r,a) behaves like some 
power of r. Inserting eq. ( P5| ) into eq. (|T2p we obtain the asymptotics of the amplitude 



at 



p ^ oo, g ^ oo, 



9= 1 



fixed, 



P 



X 



1-2k^ = fixed 



in the following form 



A, 



^(2p+l)\{2q)\ .\{yq/p+l) 



2p 2q 1 

rrii m2 



q+p 1 



- C{p, q) 

yq 



(29) 



]K'^{q/p — 1)- 

where C(p, q) behaves like a power of p and g, ?/(6', x) is given by eq. (|2^) . Obviously, 
this expression has the form of eq. (|l^ . 

Eq. (P^) looks singular at p/g — > 1, i.e. 6 = 2. However, the singularity is, in 
fact, absent because y = — 1 at that point. The asymptotics of the amplitude in this 
regime can be obtained either by the study of the integral in eq. (^) or by making use 
of eq. (|l^) and the appropriate asymptotics of the hyper geometric function. Namely, 



at A — s> oo one has 19 
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where 



and C(A) behaves as some power of A at A ^ oo. Eqs. (pOD and ([T^) lead to eq. (|T5D 
in a straightforward way. 

In this Appendix we did not discuss the pre-exponential factors in equations hke 
eq. p^) or eq. (|30|). In principle, there could have occurred cancellations between dif- 
ferent terms in eq. (0) or (p^S]). We have checked by calculating the pre-exponentials 
that the cancellations occur in none of these formulas, i.e., the exponential behaviour, 
eqs. (0), ( 1T5|) indeed takes place. 
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Fig. la 



Fig. lb 



Figure 1: x — > 3% diagrams in the model with broken reflection symmetry (/? 
solid and dashed lines correspond to the fields % and 99, respectively. 






< 




Fig.2 



Figure 2: Diagrams for the decay </? ^ 4% in the model with broken reflection sym- 
metry. 
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